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(Dated:) 

We study thermoelectric transport at low temperatures in correlated Kondo insulators, motivated 
by the recent observation of a high thermoelectric figure of merit(ZT) in FeSb2 at T ~ 10A^. Even 
at room temperature, correlations have the potential to lead to high ZT, as in YbAl^, one of the most 
widely used thermoelectric metals. At low temperature correlation effects are especially worthy of 
study because fixed band structures are unlikely to give rise to the very small energy gaps E g ~ 5kT 
necessary for a weakly correlated material to function efficiently at low temperature. We explore 
the possibility of improving the thermoelectric properties of correlated Kondo insulators through 
tuning of crystal field and spin-orbit coupling and present a framework to design more efficient 
low-temperature thermoelectrics based on our results. 



I. INTRODUCTION 

Thermoelectrics support a voltage drop in response 
to a modest temperature gradient. Since a temperature 
gradient affects the electrons and the lattice degrees of 
freedom, optimizing thermoelectrics involves not only the 
thermopower or Seebeck coefficient (S), but also the elec- 
trical (cr) and thermal (k) conductivities. The holy grail 
of thermoelectrics is to achieve a figure of merit 

zr=(^)r (i) 

that exceeds unity at room temperature. This tall order 
remains a grand challenge problerrPP Two promising 
recent directions have focused on either decreasing the 
thermal conductivity as in the case of nanocrystalline ar- 
rays of Bi a; Sb2- 3 ;Te3 in which a ZT of 1.4 was achievecP 
at T — 373K or maximizing the numerator of Z through 
strong electron correlations. An example of the latter is 
the reporlP that FeSb2 achieves a colossal thermopower 
of 45000^y/ii at 10K resulting in the largest power fac- 
tor, S 2 a witnessed to date. In this paper, we follow-up 
on the role strong correlations play in maximizing the 
power factor by focusing on Kondo insulators. We show 
explicitly that multi-orbital physics in Kondo insulators 
lies at the heart of the problem of maximizing the power 
factor. 

Because the thermopower is related to the entropy 
per carrier, particle-hole asymmetry and large density of 
states at the chemical potential are central to the op- 
timization of Z. In this regard, the Anderson model 
of a single impurity in a metaP, which is among the 
few solvable strongly correlated systems solvable exactly, 
presents a density of states with demanding features for 
efficient thermoelectric transport. For a single SU(2) 
spin on a localized impurity, the density of states appears 
as a single infinite symmetric peak at the chemical poten- 
tial leading to a divergent density of states but vanishing 
Seebeck coefficient by virtue of the particle-hole symme- 



try. Increasing the degeneracy of the localized orbital and 
the metallic band to SU (N) (N > 2) softens the peak in 
the density of states and at the same time moves above 
the chemical potential leading to an asymmetric density 
of states and as a result a larger Seebeck coefficient!-^. 

It makes sense then to consider systems in which such 
physics is naturally present, for example Kondo insula- 
tors in which a regular lattice of Anderson impurities 
is hybridized with multiple bands of itinerant electrons. 
The electrons in the local orbitals are poorly screened 
and the strong Coloumb repulsion prohibits them from 
being multiply occupied. Contrary to the single impurity, 
the periodic Anderson model is not exactly solvable but 
multiple mean-field type methods were usecP^ i un- 
derstand many of their features. Motivated by the single 
impurity model, we examine the effect of degeneracy of 
the local impurities and the conduction band on the ther- 
moelectric properties of Kondo insulators. In addition to 
directly studying the degeneracy of the local and conduc- 
tion bands, we study the effect lifting the degeneracy by a 
crystal-field (which mainly effects the local orbitals) and 
spin-orbit coupling (which mainly effects the conduction 
band) have on thermoelectric efficiency. In this way, we 
can continuously lift the level of degeneracy. Interest- 
ingly, we observe that there is an optimum value of the 
crystal field and spin orbit coupling. As was shown in a 
previous study^, the presence of multiple orbitals close 
to the chemical potential is a common feature of Kondo 
insulators. Our results indeed present a possible route for 
using strong correlations to enhance the thermoelectric 
performance through controlling the orbital degeneracy 
of local and itinerant bands. 



II. MODEL AND METHODOLOGY 

Heavy fermion materials typically contain rare earth 
or actinide ions forming a lattice of localized magnetic 
moment^l. The strong Coulomb repulsion of electrons 
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localized in / or d orbitals leads to formation of these 
local moment^ which then hybridize with the itinerant 
electron bands and form the heavy electron bands. If 
the chemical potential is in the heavy electron bands, a 
heavy fermion metal is formed. The volume of the Fermi 
surface in this correlated state corresponds to a sum of 
the number of itinerant and localized electrons. If the 
chemical potential is in the hybridization gap, the heavy 
electron ban d wil l be fully occupied and a Kondo insu- 
lator obtaina^EJ] Notice that such an insulating state 
is fundamentally different from a non-interacting insula- 
tor. For example, in order to reach a filled valance band, 
we need to add the number of localized and itinerant 
electrons which is solely developed as a result of strong 
interactions. 

The underlying microscopic model of this correlated 
system is: 
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where cL (dL ) is the creation of a conduction(local) 
electron with momentum k, orbital /, and spin a = (\, \,), 
and nf la — d., du^ is the number operator of a local 
d orbital at site rj. The dispersion of the c-electron 
e aa > (k) = £k<5o-(T' + Tk • cr aa i includes spin coupling. The 
non-dispersive energy of local states (e/,) depends on the 
orbital index I. The pseudovector represents the am- 
plitude of the spin-orbit (SO) coupling^ED an( j its form 
depends on the crystal symmetry of the underlying lat- 
tice (see appendix[A| . Typically, the hybridization matrix 
element, V\aa encodes the complex orbital structures of 
local states which can have novel effects on the properties 
of the strongly correlated heavy fermion phase^, but as 
in other studies, we consider V\a a to be independent of 
(k, a) to make the calculation more tractable. 

Using the model Hamiltonian, Eq. |2]), we can cap- 
ture the effect of the degeneracy of both localized and 
itinerant bands, as well as the effect of crystal field and 
spin orbit coupling in breaking the degeneracy of these 
bands. As a result of weak screening of electrons in / 
and d orbitals, the associated on-site repulsive potential 
U is much larger than the hopping energies of the itiner- 
ant electrons. To treat the large on-site repulsi on t erm, 
we use the U(l) slave-boson mean- field theory^E^. In 
this treatment, the creation operator of a local electron 
d\, = fjio-bi is partitioned into a neutral fermion fj l(r , 
and a charged boson hi that accounts for annihilation of 
an empty state. Since the local Hilbert space is restricted 
to either an empty or a singly occupied state, the addi- 
tional local constraint, 
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should be enforced at every site r^. The Hamiltonian in 
terms of these slave particles then becomes 
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where A^ is a Lagrange multiplier to maintain the local 
constraint. In the above Hamiltonian, the effect of the 
crystal field is to break the degeneracy of the local or- 
bital states whereas the spin-orbit coupling breaks 
the spin degeneracy of the conduction band. As a re- 
sult, by tuning the crystal field and spin-orbit coupling, 
we can change the degeneracy of the local and conduc- 
tion orbitals in a continuous manner. Consequently, we 
have a tunable knob to gain the optimum thermoelectric 
performance. 

The mean field approximation to the model Hamilto- 
nian can be obtained by taking the coherent expectation 
b = (bi) = (b\) and A = (A;). This replacement effectively 
renormalizes the mixing matrix element V\ — > bVi, and 
the local energy e/ ; — > e/, + A and leads to the quadratic 
Hamiltonian 
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Instead of working in the spin basis, we will use a he- 
lical basis that diagonalizes the single-electron disper- 
sion £ CTCT /(k) -)• Ule(k)U k = (e k + h\T k \)5 hh > with 

L J hh' 

h,h' = ±1. Then c\ah{fkih) is accordingly rotated by 
the unitary matrix from the spin basis, Cki CT (/kio-)- 
By performing the Bogoliubov transformation, 
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we obtain the diagonal mean-field Hamiltonian, 

H MF = T, E± h al h+ au h ± + A]T(& 2 - 1), (8) 
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where the dispersion is given by 
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Minimization of the free energy with respect to the 
mean field parameters b and A, and the total chemical 
potential [i leads to two coupled equations 

1 = b 2 + al lh n F {E Mh+ ) + pl lh n F {E uh _) , (12) 
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where, the total density of electrons is fixed to be ntot = 
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The transport of this non- 



interacting mean-field Hamiltonian is now tractable. 

To compute the transport properties, we use 
the relaxation-time approximation to the Boltzmann 
equation^. Under this scheme, the electrical resistiv- 
ity, p — cr^ 1 , and the thermopower tensors, S, are given 
by 
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explicitly. Here we set ~v\ah± = <fV-Ekifc±. Consid- 
ering that the electrons are scattered by -/V imp impuri- 
ties with an interaction strength of Vj mp (i.e., H sct t ~ 
Ckicr), the relaxation time T\ah± for each state is 
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with pih{E\sih±) the density of the states of the Bogoli- 
ubov quasiparticles. 

III. RESULTS 

We now present our results on the dependence of the 
transport properties on the orbital degeneracies of both 
localized and itinerant electron bands which form corre- 
lated Kondo insulators. In the first two subsections, we 
consider double degeneracy of conduction and localized 
bands. This model is indeed consistent with the models 
previously proposed for Kondo insulators 2 ^. The crystal 
field will then split the degeneracy of the two / levels 
in [4] and spin-orbit coupling breaks the degeneracy of 
the conduction band states with differing helicitjffS In 
these two sections, we change the size of the degeneracy- 
breaking gap continuously. Using the relaxation-time ap- 
proximation, we can then calculate the transport prop- 
erties of a Kondo insulator. For most of the materials, 



the dominant contribution to the thermal conductivity 
comes from lattice vibrations; as a consequence, the elec- 
tronic contribution to the thermoelectric performance is 
measured through the power factor Zpp = aS 2 where a 
is the electrical conductivity and S is the Seebeck coef- 
ficient. In order to confirm that the enhancement of the 
thermoelectric efficiency can properly be attributed to 
strong correlations, we consider two different band struc- 
tures of itinerant electrons in IHI Al and IIIIBI and we see 
that similar features emerge. 

Finally, in |III C| we present the effect of multiple orbital 
degeneracy. Contrary to the treatment in |III A| and |IIlB| 
where the double degeneracy is continuously lifted by 
crystal field and spin-orbit couplings, in section [Til C| we 
discretely change the number of degenerate conduction 
and localized bands. We show that indeed there is also 
an optimum degeneracy associated with the maximum 
power factor. 



A. Nearly free electron itinerant bands 

We first focus on the effect of crystal field and spin- 
orbit coupling on the power factor within the context 
of a parabolic band for the itinerant electrons ek = 
e + W(k/k BZ ) 2 with W = 2eV, taken from refP. 
In principle, the SO coupling should be expressed as 
a periodic function under the crystal environment, but 
we model it to be isotropic as well, Tk = 7so(k/^sz) 
(7so < 0.2 eV). In the following, we carry out the nu- 
merical calculation based on this isotropic band disper- 
sion with Z max = 2. Here, we choose ej 1 = 1.0606eV, 
Vi = 0.2236eV and V 2 = 1.05Vi = 0.2348eV. The pa- 
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0.045eV 5 EI]. The other control 
parameters are temperature(T < 100AT), crystal electric 
field (CEF) splitting(A C EF = e/ 2 - e fl < 15meV), and 
the SO coupling(7 so < 0.15 eV). Although we do not 
include the supporting data here, we found that our con- 
clusions are insensitive to the strength of as long as 
0.5 < V 2 /Vi < 2.0. 

Figure [l] shows the density of states (DOS) for dif- 
ferent control parameters. From Figs, [lja) to (d), we 
notice that the temperature only controls the number of 
thermally activated charge carriers, while it does not sig- 
nificantly change the DOS compared to the other param- 
eters. When the degeneracy of the two local orbitals is 
broken by the crystal electric field, one of the hybridized 
bands moves closer to the chemical potential. Conse- 
quently, the system is driven from an insulator to a con- 
ductor [Figs, [lje)-(h)], at which point the power factor 
is significantly enhanced (See Fig. [2]). Likewise, the spin- 
orbit interaction breaks the degeneracy of the two heli- 
cal modes, which turns an insulator into a metallic state 
[Figs, nji)-(l)]. We point out that the metallic state is 
characterized either by a local orbital, I = 2 [Fig. |T|(h)] 
or by a helicity, h = + [Fig. [ljl)], since only the bands 
with corresponding quantum numbers are conducting. 
In order to further examine the effect of the control 
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FIG. 1. Density of states (DOS) of the parabolic model for 
different control parameters. The red-dashed lines are for the 
orbital 1 = 1, the dotted blue lines I = 2, and the green 
lines are for the total DOS. The central gray area displays 
the thermal window(~ ksT) for each temperature. Figures 
(a)-(d) compare the DOS for different temperatures. Figures 
(e)-(h) correspond to different magnitudes of the crystal elec- 
tric field(AcBF) which breaks the degeneracy of the two / 
orbitals. Figures (i)-(l) correspond to different magnitudes of 
the spin-orbit interaction(7so) that breaks the degeneracy of 
the conduction bands with different helicity. 
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FIG. 2. Transport coefficients for different temperatures as 
a function of crystal field splitting and spin-orbit interaction: 
thermopower S (left column), conductivity a (middle), and 
power factor Zp? (right). From the top to the bottom, the 
temperature varies from 50 to 5K, and the solid lines are 
equally spaced constant contours. 



parameters, we first calculate the transport coefficients 
as a function of the CEF and SO. In Fig. [2] we show 
the results of a calculation of the thermopower (S) , the 
electrical conductivity (a), and the power factor (Z PF ). 
As can be seen from the right column, the power fac- 
tor is enhanced either by finding the optimal CEF or by 
adjusting the SOI. Since both CEF and SOI shift some 
of the lower energy bands toward the chemical potential 
[Figs. [TJe)-(l)] , the number of lower energy bands rele- 
vant for thermal transport is controlled by CEF and SOI 
simultaneously. For the temperature range T < 20K, 
where the thermal windows are sufficiently narrow, CEF 
and SOI compete; hence there are two distinctive op- 
timal regimes. For a sufficiently wide thermal window, 
attainable at intermediate temperatures, T ~ 30-fT, CEF 
and SOI are working cooperatively to form a single op- 
timal region. For T > 30K, the enhancement in Zpp is 
not as drastic as at low temperature. Zpp is maximized 
in the vicinity of the insulator-metal transition (see the 
conductivity a at T = 5 — 20-ftT), resulting from a com- 
petition between S and a. For instance, at T = bK, 
the thermopower S decreases with SOI and CEF, while 
the system acquires a finite conductivity. Note that the 
metallic state here has one dominant helical state over 
the other. 

In Fig. [3j we repeat the calculation of the transport 
coefficients for a fixed SOI as a function of CEF and 
temperature. Consistent with Fig. [2] is that the optimal 



point for the power factor is located in the vicinity where 
the insulator-metal transition occurs. For instance, when 
7so < O.leV, there is the optimal CEF and temperature 
for the power factor, at which point the electric conduc- 
tivity acquires a noticeable finite value. From the left 
column, one finds that the thermopower generally de- 
creases with increasing temperature as a widened ther- 
mal widow implies the reduction of the asymmetry in 
the DOS within the thermal region [Figs. [TJa)-(d)]. This 
obtains because as the temperature increases, more of 
the bands (lower and upper) are involved in the thermal 
transport. In other words, the asymmetry of the DOS 
within the thermally active region is relieved. Beyond a 
certain of the threshold SOI, 7 so > 0.13eV, there is no 
phase transition (at mean field level) ; hence optimization 
cannot be realized. 



B. Tight binding itinerant electron bands 

Next, we consider the 3-dimensional tight binding case. 
We see that as in the case of a quadratic band, tuning 
the crystal-field and spin-orbit coupling can optimize the 
thermoelectric performance. This result indicates that 
the effect of orbital degeneracy in controlling the ther- 
moelectric performance is not that sensitive to the de- 
tails of the band structure. Here, we choose the hopping 
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model: the power factor can be enhanced by adjusting 
the CEF, while the SOI slightly lowers the optimal value 
of the CEF. For T > 15K, the enhancement in Zpp is 
not as drastic as was in the low temperature case. As in 
the parabolic model, this trend occurs because the ther- 
mally active region is too wide to encompass only one 
band [see Fig. [4^c) and (d)]. Comparison with the other 
columns reveals that Zpp is also maximized near an insu- 
lator to a (helical) metal transition (see the conductivity 
a at T = 5K), which is the consequence of the compe- 
tition between S and a. Here, one can observe that S 
becomes maximal at Acef — O.Olihop, which is a conse- 
quence of the choice V 2 /V 1 = 1.01. With V 2 /Vi = 1, S 
only decreases with CEF (not shown) . Fig. [6] similarly 
confirms the consistency with the parabolic model. The 
only difference is that the effect of SOI is not as remark- 
able, though it works to shift the optimal value of the 
CEF. The reason mainly lies in the changes of the DOS 
depending on SOI: linearized SOI changes the bandwidth 
significantly, while 3D tight-binding SOI does not due to 
its periodic structure. (See Fig. [7]) 



FIG. 3. Transport coefficients for each SOI as a function of 
CEF and temperature: thermopower S (left column), conduc- 
tivity a (middle), and power factor Zpp (right). The range of 
SOI is 0.02 — 0.15 eV from the top to bottom panels. 



parameter t^ op = 0.2167eV (band width W = 2.6eV), 
and we located the local energy €f 1 = — 0.8iho P - The hy- 
bridization strength V\ = ihop and U 2 — l.Olthop- In the 
SOI, we take the next-nearest neighbor hopping param- 
eter gi = 0.3. 

As in the simplified parabolic model, the role of CEF 
and SOI is not different; both efficiently control the sys- 
tem to drive it from an insulator to a conductor as seen 
from Fig. |4j Compared to the corresponding panels in 
Figs, [l] however, Figs. llFe)-(h) show that the CEF also 
pushes one of the upper bands toward the chemical po- 
tential, hence reducing the gap size significantly. In fact, 
the parabolic model is rather exceptional since the bot- 
tom of the upper bands corresponds to the point k = 0, 
which is not usual for typical 3D tight-binding models. 
Figs. [2ri)-(l) display the evolution of DOS with the in- 
crease of the SOI. Even though the degeneracy of the two 
helical modes are broken with a finite SOI, it cannot be 
seen clearly as was in the linearized SOI case [Figs. [TJi)- 
(1))]. The reason is that \Tk\ decreases as k approaches 
the boundary of the Brillouin zone due to the periodic 
form of the SOI, while it does not for the linearized SOI. 
Unlike the CEF which affected drastically only one of 
the orbitals, the effect of the SOI is quite different. Up 
to 7so = 0.2i hop ~ 40meV, the changes in the DOS are 
not significant. For 7 so > 0.2ihop, the system undergoes a 
phase transition to a (helically polarized) metal, beyond 
which point Zpp is reduced (See Fig. [5]). 

From Fig. [5] we observe consistency with the parabolic 



C. Effect of multi-orbital degeneracy 

In addition to the continuous control of orbital degen- 
eracy through crystal-field and spin-orbit coupling, we 
can specifically study the effect of increasing the number 
of degenerate orbitals (/ max = 1, 2, • • • ,5). To minimize 
the number of free parameters, we will set the orbital de- 
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FIG. 4. Density of states(DOS) of the 3D tight-binding model 
for different control parameters. Both Acef and 7 so are in 
units of the hoping amplitude, t^ op = 0.216eV. The red- 
dashed lines are for the orbital 1 = 1, the dotted blue lines are 
for I — 2, and the green lines are for the total DOS. The cen- 
tral gray area indicates the thermally active region for each 
temperature. Figs, (a)-(d) compare the DOS for different 
temperature, (e)-(h) for the crystal electric field(AcEF), and 
(i)-(l) for the spin-orbit interaction(7so)- 



6 




0.00 



-V 1 1 








i 









£ 

□ 

o.5 g 




0.00 0.04 0.00 0.04 0.00 0.04 

AcEF(thop) 

FIG. 5. Transport coefficients for each temperature: ther- 
mopower S (ieft coiumn), conductivity a (middie), and power 
factor Zpf (right). From the top to the bottom, the temper- 
ature is fixed to 25, 15, 10, and 5K, respectively. 
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FIG. 6. Transport properties are compared for each SOI: S 
(left column), conductivity a (middle), and power factor Z-p-p 
(right). The spin-orbit couplings are chosen to be 0.0, 0.1, 0.2, 
and 0.5 in units of the hoping amplitude thop, respectively. 



generacy of the two bands (the allowed values of /) to be 
equal. Although continuous control is not possible in this 
case, one can then consider changing the material con- 
tent to achieve a better thermoelectric. Here, the bare 
conduction electron dispersion is taken to be that of the 
3D tight binding model. 

First, we compare the DOS depending on the number 
of orbitals involved [Fig. [8]. As Z max increases, the asym- 
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FIG. 7. Density of states for bare conduction electrons: 
quadratic dispersion(top), 3-D tight binding(bottom), helic- 
ity h = — l(left, blue curves) and h = l(right, red lines). For 
the quadratic dispersion, the SOI is taken to be linear in mo- 
mentum, as 7 so |k|. The changes in the bandwidth are exactly 
proportional to 7 so . In figures (c) and (d), the SO term, 7 so rk 
is taken in accordance with the cubic point group symmetry, 
and the next nearest neighbor hopping parameter gi — 0.3. 
Note that the bandwidths are not drastically affected by the 
SO, while the shapes become more asymmetric with the in- 
crease in 7 so . 



metry between the upper and the lower bands becomes 
more pronounced. At the same time, the insulating gap 
increases with Z max for Z max > 2. Note that this feature 
is quite similar to the single impurity problem with N 
flavors. The inset of Fig. [8] displays the DOS without 
adjusting the chemical potential. Since the slave-boson 
method renormalizes the local energy by e/ — > e/ + A, 
the relative location of the Kondo resonance for each 
case (near each gap) indicates that the amount of renor- 
malization A increases with the number of available or- 
bitals. Since each band below and above the insulating 
gap should accommodate one electron, the deformation 
of the lower bands becomes less significant as ^ max in- 
creases (see the lower bands for different / max ). In other 
words, since the area below and above the gap should be 
equal, the asymmetry of the DOS becomes more signifi- 
cant as A, or equivalently l max , increases. 

Given the band structures at the mean-field level, we 
proceed to calculate the transport properties as shown 
in Fig. [9] As / max increases, maxima of Zpp also in- 
creases as the temperature is elevated. The thermopower 
is also enhanced with the number of available orbitals, 
which is caused by pronounced asymmetry in the DOS 
(See Fig. |8j). Since the gap size increases, the conduc- 
tivity generally decreases with the number of orbitals 
for i max > 2. Though not shown here, the maximal 
power factor per orbital, Zp-pJJm&x, also increases with 
'max until Z max < 7. In Fig. K5lc), the thermal conduc- 
tivity due to electronic structure is evaluated, exclud- 
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FIG. 8. Density of states per orbital for different number of 
orbitals(£ max ), where the chemical potential at each case is 
adjusted to the center. The inset displays the DOS without 
the adjustment. 
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FIG. 9. Transport properties for Z max orbitals: (a) Seebeck 
coefficient, (b) electrical conductivity, (c) thermal conductiv- 
ity, (d) thermopower (e) figure of merit and (f ) dimensionless 
figure of merit. 



ing any contribution from lattice vibrations. Typically, 
phonons are dominant contributors to the thermal con- 
ductivities, but it may not be so prevalent at the low tem- 
perature range considered here, presumably T < 10K. 
The resultant (dimensionless) figure of merit, assuming 
K = ^electron + K phonon — n — ^electron? is strongly en- 
hanced with the number of orbitals at least by an order 
of 10. This is one of the key results of this paper. 

IV. SUMMARY 

We have shown that orbital degeneracy opens a new 
window into increasing the power factor for strongly in- 
teracting thermoelectrics. Our key findings are that the 



power factor is maximized in strongly correlated systems 
by tuning 1) the gap between nearly degenerate local 
f-orbitals through the crystal field effect, 2) ) the gap be- 
tween nearly degenerate itinerant electron bands through 
the spin-orbit coupling, and 3) the the number of degen- 
erate local and itinerant orbitals. 

This approach provides a new parameter space for the 
design of strongly correlated thermoelectric materials. 
Our result was derived using the slave-particle mean-field 
theory, which is not expected to be quantitatively reliable 
but should capture general trends. The effect of degener- 
acy in enhancing thermoelectric performance of strongly 
correlated systems could also be investigated using other 
methods such dynamical mean-field theory^ and finite 
frequency methods^. Another direction for future work 
is to consider more profound effects of spin-orbit coupling 
combined with correlations, as in the proposed "topolog- 
ical Kondo insulators'^, whose surface states are cur- 
rently being sought experimentally; such surface states 
have the potential to increase thermoelectric performance 
at low temperatures^. 

Finally, a direction that is difficult theoretically but 
may be important for actual materials is to find ways of 
interpolating between the effectively itinerant calculation 
here (i.e., there are plane- wave states of the slave parti- 
cles) with the "atomic limit"^, where the effects of mul- 
tiple orbitals have also been considered^. The atomic 
limit, which is valid when the hopping is the smallest 
energy scale in the problem, has been argued to be rel- 
evant to ex perim ents on sodium cobaltates near room 
temperatur e 1 27 ! 28 !. The results of this paper should mo- 
tivate continued investigation of correlated materials for 
thermoelectricity and suggest that a guided search with 
controlled crystal-field splitting may lead to further im- 
provements in thermoelectric figure of merit, especially 
in the low-temperature regime. 
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Appendix A: Spin-orbit coupling 

In the presence of SO coupling} 15 * 16 !, fa e conduction 
electron dispersion matrix e aa i (k) is given by 

£ CTCr < (k) = e k <W + T k • tr aa > , (Al) 

where ek is the dispersion without the SO interaction, 
and (7 are the Pauli matrices. (The SO interactions con- 
sidered here originate from the absence of an inversion 
symmetry in the crystal lattice.) The antisymmetric SO 
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coupling is described by the real pseudovector Tk, which 
is determined by the point group symmetry of the crys- 
tal. For instance, CePt 3 Si, CeRhSi 3 and CeIrSi 3 belong 
to tetragonal point group (G = C± v ) in which 



Tk = 7so(kz sin k y a - k y sin k x a 



,(A2) 



+ k z <7 2 sin k x a sin k y a sin k z c(cos k x a — cos k y a) 

in the next nearest neighbor approximation for a real 
7so, the lattice spacing a, c and the next nearest neighbor 
parameter g 2 - In case of the cubic point group symmetry 



(G = O), the pseudovector is given by 

T k = -f so k x smk x a[l - g 2 (cos k y a + cos k z a)] 
+ (Positive permutations of x, y, z). 

In real noncentrosymmetric crystals, the typical SO 
strength ranges up to 200 meV. Instead of working 
in the spin basis, it is useful to introduce the helical 
basis that diagonalizcs the single-electron dispersion 



<w(k) -5 
h,h' = ±1. 



Ule(k)U k 



hh> 



(e k + h\T k \) 5 hh > with 
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